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Radiative emission from excited states of atoms and molecules can comprise a 
significant fraction of the total heat flux experienced by spacecraft during 
atmospheric entry at hypersonic speeds. For spacecraft with ablating heat shields, 
some of this radiative flux can be absorbed by molecular constituents in the 
boundary layer that are formed by the ablation process. Ab initio quantum 
mechanical calculations are carried out to predict the strengths of these emission 
and absorption processes. This talk will describe the methods used in these 
calculations using, as examples, the 4 th positive emission bands of CO and the 1 2 g + 
1 2 u + absorption in C3. The results of these calculations are being used as input to 
NASA radiation modeling codes like NeqAir, HARA and HyperRad. 
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Introduction 
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Modeling radiative heat flux fo hypersonic 
flight (~10 km/s, 50-90 km altitude) 

- Spectroscopic databases needed for NASA's 
existing radiation codes (NEQAIR, HARA) and 
new codes under development (e.g., 

HyperRad) 

- These codes use computed flow fields to 
predict radiative transport in the shock and 
boundary layer during re-entry 

- Validation using ground facilities (EAST, 

LENSxx, X3) 

- Ultimate goal: accurate prediction of heat flux 
on space vehicles during re-entry in order to 
reduce safety margins on thermal protection 
system 




Re-Entry Conditions 
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Earth (N 2 /0 2 ) 

- Hypervelocity return (> 9 km/s) 

- Large blunt bodies (several meter diameter) 

- Ablating heat shields (carbonaceous species injected 
into boundary layer) 

Mars (C0 2 /N 2 ) 

- ~6-8 km/s 

- MSL-size bodies and larger 

- Aerobraking at high altitudes (large flexible heat shield) 

- Conventional heat shields at lower altitude and speed 


Important radiative processes 
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Most radiation is from atomic emission (neutrals 
and cations) 

Strongest molecular emission from VUV bands 

- Much of this is re-absorbed by shock and boundary 
layer species before it reaches spacecraft surface 

Also contributions from bound-free and free-free 
continuum radiation (not within the scope of this 
presentation) 


Molecular Species of Interest 
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Earth 

- N 2 , 0 2 , NO, N 2 + , NO + (emission in shock layer) 

- C 2 , CO, C 2 H, C 3 H, C 3 , C 4 , etc. (absorption in boundary 
layer) 

Mars 

- C0 2 (IR), CO (IR, VUV), CN (UV-Vis), N 2 (UV) 


Spectroscopy: Top down or Bottom up? 
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The Old Way The New Way 



Start with experimental • 
spectrum • 

Determine initial and final 
ro-vibrational levels 

- Compute RKR potentials 

- Franck Condon factors 

Measure absolute (or 
relative) intensities of 
individual lines 

- Compute R-centroids 

Compute Einstein A- 
coefficients 


Quantum mechanics all the way! 

Electronic Schrodinger equation 

- Compute potential energy curves 
(surfaces) 

- Compute transition moments and 
coupling terms 

- Obtain analytic fits 

Nuclear Schrodinger equation 

- Compute ro-vibrational levels 

Direct determination of Einstein 
A-coefficients 


Quantum Mechanics Refresher 




Electronic Schrodinger equation (clamped nuclei, 
shown for diatomic molecules) 


V (potential) 

* electron-electron and electron-nucleus Coulomb 
interaction for potential energy calculations 

* diatomic potentials or triatomic PES for ro-vibrational 
energy level calculations 

Wavefunctions (u) and energy (E) determined by 
basis set expansion methods 



Spectroscopic intensity factors 
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Franck-Condon factors and R-centroids 




/ 


r dr 


/! 


dr 


Einstein A coefficient for spontaneous emission 

A^= (2.026 x 10' s ) % viv [JvtfKc(r)MrJ 

A^^is the probability per unit time that a 
molecule in the upper state will decay to the 
lower state by emitting a photon 


Our Objective 
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Use quantum mechanical calculations to simulate 
molecular spectra from the IR to the VUV 

- Diatomic and triatomic molecules 

- Minimize empirical adjustments 

- Complete set of initial and final ro-vibrational levels 

- Electric and magnetic dipole and quadrapole transitions 

- Weak and strong band systems 

Simulate absorption and emission spectra for 
high-speed re-entry conditions 


New Wrinkles in the Present Work 
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Diatomic potentials accurate all the way to the 
dissociation limit (not true for RKR and Morse- 
type potentials) 

Potentials include proper interactions due to 
coupling with other electronic states (e.g., spin 
orbit) 

Diabatic representation (instead of adiabatic) 
removes 'avoided crossings' by allowing states to 
cross 

Improved treatment for Rydberg states (cause of 
many of the avoided crossings) 


Avoided crossings 
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Adiabatic representation 

- diatomic potentials and PESs of 
same multiplicity and symmetry 
do not intersect as the red and 
green curves (top) illustrate 

- Avoided crossing results in 
change of electronic structure 
as shown by oscillator strengths 
(bottom) 


C 3 6 1 A', 7 1 A' & 8 L A' states 
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Diabatic representation 
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Adiabatic representation 

- Electronic states of same symmetry do not cross 

- But electronic structure of a state changes abruptly at 
geometries where two states are nearly degenerate 
(avoided crossing) 

Diabatic representation 

- Electronic structure of each state changes gradually as 
the geometry changes 


Adiabatic-diabatic transformation 
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• Transformation effected by nuclear kinetic energy 
operator 

- Two-state example: 


Tn = Y. Y, with P Ac , = — iV.ia = -i-P 


.4 a=x,y,z 


2M a 


ORa* 
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Adiabatic vs diabatic 
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Adiabatic 

- Easier to compute 

- Potentials are more accurate, but harder to use in 
calculations of ro-vibrational energy levels and spectra 

Diabatic 

- Harder to compute 

- Easier to use in subsequent calculations of molecular 
spectra 

- Requires more higher adiabatic states to have 
adequate convergence 

- Works well around potential minima, but not at the 
asymptotes 


Adiabatic, Diabatic Example 
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NO 2 n electronic 
states 

- Avoided crossings 
removed in diabatic 
representation 

- N + O asymptote not 
well described 


NO 2Pi 



■129.8 


Difficulties computing VUV spectra 
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Typically there are many electronic states lying 6- 
12 eV above the ground state 

- This results in many avoided crossings (for states of the 
same symmetry and multiplicity) 

- Transitions with large oscillator strength are often 
spread among several upper electronic states (each one 
important for a small range or bond lengths or angles) 

- Hard to map out all of this and to fit these potentials 
and transition moments 



Difficulties for triatomic molecule calculations 
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• Triatomic molecules have 3 geometric 
parameters 

- 2 bond lengths and 1 bond angle 

- Many more energy calculations needed to represent 
3-d potential hypersurface 

• Linear triatomic molecules like C0 2 and C 3 have 
lower symmetry as they are distorted 

- D ooh for 180° bond angle and equal bond lengths (£ + , 

i u + , n g , n u , etc.) 

- C s for < 180° bond angle and unequal bond lengths 
(A', A") 



Examples today 
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CO 4 th -positive emission 

- 

- 0-0 transition at 154.4 nm 

- Typically observed in 145-195 nm wavelength range 

- Shorter wavelengths obscured by black body limit 

- Largest molecular contributor to radiative flux for Mars 
entry 

C 3 VUV absorption 

- % + x ^ 

- 0-0 transition between 160-170 nm 

- Overlaps very strong N(l) lines at 174.27 and 174.52 nm 


CO 4+: Computed electronic state potentials 
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Transition between CO A 
IT and ground state 
(X X S + ) [lines with 
symbols] 

- Four triplet states perturb 
the A state potential 

- Three Rydberg states 
shown [dashed lines] 

0-0 transition at 154.4 
nm (64,748.48 cm 1 ) 
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CO 4+ 
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This transition accounts for more than 60% of the 
integrated radiative flux for typical Mars entry at 
8 km/s 
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CO 4+ Comparison 
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Measured in EAST facility at NASA Ames (Electric 
Arc Shock Tube) - pure C0 2 at 6.3 - 8 km/s 

Simulation by NEQAIR (vintage code with current 
top-down spectral database) and HyperRad (new 
code with new bottom-up database) 

Measured instrumental width applied with 
natural line broadening factors 


A n potential curve shifted up by 0.013 eV 
(105 cm 1 ) to match experimental T e 


Comparison of Spectral Databases 
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Red curved line is black body cutoff for 6085 K 
Shock speed is 6.3 km/s 


0-0 transition marked by arrow 
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CO 4+ wrap up 
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State-of-the-art quantum mechanical calculations can 
predict ro-vibrational energies to within 10 cm-1 and 
electronic excitation energies to ~100 cm-1 (~0.01 eV) 

- Spectroscopic measurements still better for obtaining precise 
wavelengths of spectral lines 

- Simple adjustment of of T e values for electronic states can yield 
accurate spectra 

For the CO 4+ transition new database has very good 
match to EAST measurement 

- Old experiment-based database comparable at shorter 
wavelengths, but errors grow for longer wavelengths (higher 
ground state rovibrational levels) 


C spectroscopy 
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C 3 is the most stable carbon molecule in the vapor 
phase 

- Product of arc heated graphite and shocked acetylene 
or methane 

- Linear molecule with very small bend frequency 

Vis-UV spectral bands (Swings) from + 

transition observed in emission and absorption 

Poorly characterized VUV absorption transition 

due to 1 S U + <- 1 E g + 

Injected into the boundary layer during re-entry 
as a pyrolysis product of ablating heat shields 


Role of C 3 during re-entry 
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Boundary layer C 3 could absorb some of the incident 
radiative heat flux originating in the shock layer 

- Pair of very strong atomic nitrogen lines at 174.27 and 
174.52 nm lie within the wavelength range of this feature 

- High-resolution spectrum only observed at T < 10 K 
(Monninger et al., J Phys Chem A 2002, 106, 5779) 

- Possible feature seen in high-temperature shock tube 
experiment (Shinn AIAA-1982-1189) 

Objective - characterize C 3 VUV feature for T ~ 3000 
K and assess its 'radiation blocking' potential 


C 3 Experimental Low-T Spectra 
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Monninger et al. absorption 
spectra 

- Upper 10 K in solid Ar 

- Lower 4.3 K in solid Ne 

Peak intensity 160-170 nm 

- Band extends from ~145-180 nm 

C 3 generated from arc-heated 
graphite and trapped in 
cryogenic matrix 




Figure 4. Absorption spectra of the C< ('Z + u — X ‘Z + g ) system (a) in 
solid argon at 10 K and (b) in solid neon at 4.3 K. The spectra are 
recorded at a resolution of 1 A, and no background correction was 
applied. 


Experimental High-T Spectra 
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Shinn (NASA Langley) 1 % C2H2 or 2% CH4 in Ar in 
shock tube at T = 3600-4100 K 
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Fig. 6 Comparison of cross section measured 
at temperature 3650-3850 K range using 
1 -percent C 2 H 2 plus 99-percent argon 
mixture and 2-percent CH^ plus 
98-percent argon mixture. 
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Fig. 8 Comparison of C, absorption cross section 
between tests performed at two different 
temperatures with a mixture of 2-percent 
CII. plus 98-percent argon. 


C ground electronic state 
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Very small bending frequency (co e = 63.4 cm 1 ) 

- Large bend amplitude (± 60°) 

At cryogenic temperatures, only the lowest 
rovibrational energy level is populated 

At 3000 K (typical boundary layer 
temperature) hundreds of ro-vibrational 
levels are populated 

- Lowest level accounts for ~2% of the total 
population 


Radau radius 


Contour plot of ground state PES in terms of 
bending and symmetric stretch 
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C3 



C 3 Bending potential and energy levels 
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C-C-C angle, degrees 


C excited electronic state potentials 
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Next two slides show C 3 excited electronic states 
in terms of the bending and asymmetric stretch 
coordinates 

The VUV transition involves states that are 6-8 eV 
above the ground electronic state 

The Swings bands (vis-UV) involves the lowest 
excited state ( 1 n u which becomes 2 X A' + 1 1 A" at 
non-linear geometry) 

- Splitting due to Renner-Teller coupling (blue solid and 
dashed curves) 



Rel. Energy (eV) 


C bend: adiabatic vs diabatic 
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C3 Adiabatic bend 1A’ 



C3 Diabatic bend 1A' 
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asymmetric stretch: adiabatic vs diabatic 
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Computed and observed 4.3 K spectra 
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C 3 VUV spectrum at 3000 K 


Of 
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PRELIMINARY RESULT 


C3 with T=3000K 



C wrap up 
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• Computed cryogenic and 3000 K absorption spectra for C3 
VUV transition 

- At low temperature the computed spectrum agrees well with 
experimental result of Monninger et al 

- At high temperature the computed spectrum has the same shape as 
the result from Shinn's shock tube experiments, but the peak is red 
shifted to 172 nm 

• The population of many ground state ro-vibrational levels 
washes out structure seen in the 4.3 K spectrum 

- Peak in absorption spectrum is red shifted at high temperatures due 
greater contribution from higher ground state ro-vibrational levels 


Supports hypothesis that C 3 VUV can absorb N 174 nm 


Gallery of Computed Diatomic spectra: CO 
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CO lO.OOOK 



el = electric dipole 
e2 = electric quadrapole 
ml = magnetic dipole 


absorptivity (a.u.) absorptivity (a.u.) 


Gallery of Diatomic spectra: CN & C 2 
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CN 1 O.OOOK 
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Gallery of Diatomic spectra: N2 & NO 
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N2 1 0,000 K 



wavelength ( n m ) 


NO 1 O.OOOK 
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Summary 
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Bottom-up quantum mechanical approach can 
yield accurate VUV spectra for diatomic and 
triatomic molecules at high temperatures 

Only small empirical adjustment of upper 
electronic state T e values needed 

For CO 4+ the resulting spectra are in somewhat 
better agreement with spectra derived from EAST 

For C 3 resulting VUV spectra are much more 
accurate than anything previously available 

Results support hypothesis that C 3 VUV can 
absorb N 174 nm 


Final thoughts 
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Ab initio molecular VUV spectroscopy is competitive 
with traditional approaches for diatomics and 
represents a major advance when it comes to 
triatomics. 


Studies should be carried out to assess the 
effectiveness of C 3 in blocking N 174 nm radiation. 
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